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Abstract— Este trabalho apresenta uma nova abordagem na 
estimação de harmônicos e inter-harmônicos em sinais de cor- 
rente e tensão em sistemas elétricos de 13,8 kV e 69 Kv. 
A abordagem é alternativa às metodologias que utilizam a 
transformada rápida de Fourier (FFT). Ela combina dois es- 
timadores paramétricos: o estimador de Prony modificado e o 
estimador de Kalman. O estimador de Prony modificado fornece 
as fregiiências, enquanto o estimador de Kalman fornece a 
evolução das amplitudes dos componentes senoidais presentes 
no sinal analisado. A metodologia permite que as distorções nas 
formas de onda sejam monitoradas por um longo período. 


Index Terms— Estimação paramétrica, harmônicos, inter- 
harm”nicos. 


I. INTRODUÇÃO 


O mercado de energia elétrica vem experimentando um 
processo de desregulamentação nos últimos anos. É previsto 
que no futuro isso gere um ambiente de competição entre as 
empresas do setor. 

O avanço da tecnologia tem conduzido à instalação de car- 
gas mais sensíveis a distorções da rede elétrica, por outro lado, 
fiscalizações rigorosas e com níveis de exigência maiores são 
adotadas pelas agências reguladoras sobre as concessionárias 
de energia. Assim, é natural que as normas ou recomendações 
que tratam de questões sobre qualidade de energia se tornem 
gradativamente mais importantes para as concessionárias. 

As normas vigentes, como a [1], ou a [2] recomendam 
limites tanto para as distorções harmônicas da tensão fornecida 
pela concessionária, quanto para as injeções de correntes 
harmônicas e inter-harmônicas provenientes dos clientes que 
utilizam cargas não lineares. 

A correta estimação dos níveis de harmônicos e inter- 
harmônicos do sistema elétrico é de extrema importância 
tanto para a concessionária quanto para os seus clientes. A 
identificação de fontes de harmônicos no sistema elétrico 
depende da estimação do espectro das formas de onda de 
tensão e corrente em seus barramentos [3]. Procedimentos que 
tentam repartir a responsabilidade da concessionária e de seus 
clientes sobre as distorções harmônicas, também necessitam 
da estimação dos harmônicos [4], [5]. 

O conhecimento dos harmônicos e inter-harmônicos no 
sistema elétrico também permite que as concessionárias ela- 
borem políticas que incentivem seus clientes a limitarem a 
geração de correntes harmônicas. 6 realizaram um estudo sobre 
várias possibilidades de estruturação tarifárias que penalizam 
economicamente os consumidores que geram distorções nas 
formas de onda de corrente e de tensão. Por outro lado, os 


consumidores necessitam que o espectro das formas de onda 
de corrente geradas por seus equipamentos seja conhecido. 
Dessa forma, filtros mais eficazes podem ser projetados. 

Em relação aos inter-harmônicos, é sabido que suas 
frequências e amplitudes são variantes no tempo, dificultando 
ainda mais a análise dos sinais [7]. O desempenho de filtros 
passivos na presença de inter-harmônicos é limitado [8]. A 
atuação desses filtros é efetiva em uma faixa fixa do espectro 
dos sinais. Uma solução é a utilização de filtros ativos. Eles 
geram e injetam no sistema sinais de correntes iguais em 
amplitude e fregiiência, mas de fases opostas, a aqueles com- 
ponentes que se deseja filtrar. Sistemas de monitoramento de 
equipamentos também são baseados na estimação do espectro 
dos seus sinais. [9]. 

Tradicionalmente, a transformada rápida de Fourier (FFT 
- Fast Fourier Transform) é a ferramenta mais utilizada na 
análise de sinais da rede elétrica. Contudo, a correta aplicação 
da FFT parte das seguintes premissas [10], [11]: 

1) O sinal a ser analisado deve ser periódico. 

2) A frequência de amostragem deve ser maior que duas 

vezes a maior frequência do sinal. 

3) A freqiiência de amostragem deve ser igual ao número 
de amostras multiplicado pela frequência fundamental 
assumida previamente no algoritmo. 

4) Cada freqiiência contida no sinal deve ser um múltiplo 
inteiro da fregiiência fundamental assumida previamente 
no algoritmo. 

As premissas citadas acima podem ser traduzidas da 
seguinte maneira: Para que a FFT possa ser aplicada correta- 
mente é necessário que o produto entre o número de amostras 
(pontos) e o período de amostragem corresponda a um número 
inteiro de ciclos do sinal. Além disso, o teorema de Nyquist 
deve ser respeitado. 

Em relação aos sinais oriundos da rede elétrica, duas dificul- 
dades surgem na aplicação da transformada rápida de Fourier. 
Primeiramente, as amplitudes dos componentes harmônicos e 
inter-harmônicos presentes no sinal analisado podem ser vari- 
antes no tempo. Assim, o sinal analisado não seria estacionário 
e a FFT não poderia ser aplicada corretamente [12]. Além 
disso, o período fundamental do sinal nem sempre corresponde 
ao componente fundamental de 60 Hz, e a obtenção de uma 
janela de tempo que corresponda a um número inteiro de ciclos 
não pode ser realizada. 

Esta trabalho propõem um novo procedimento de análise 
espectral de sinais de tensão e corrente em sistemas elétricos. 
A análise deve fornecer a evolução das amplitudes dos 


harmônicos e inter-harmônicos presentes nos sinais analisados. 
Neste contexto, o termo evolução significa rastreamento das 
variações das amplitudes das senóides no decorrer do tempo. 
Para tanto, os sinais são coletados através de transformadores 
de corrente e de potência instalados nas subestações de entrada 
de plantas industriais ou em subestações de proteção e medição 
da concessionária elétrica. 

As amplitudes harmônicas e inter-harmônicas serão esti- 
madas através de um filtro de Kalman que é um estimador 
linear capaz de rastrear de maneira eficiente mudanças nos 
parâmetros monitorados. No entanto, é necessário que as 
fregiiências contidas no sinal analisado sejam previamente 
estimadas. 

A estimação de freqiiências é um problema complexo, 
pois é uma otimização extremamente não linear. Nesta Tese, 
desenvolveu-se uma nova técnica de estimação de freqiiências 
baseada no método de Prony e denominada método de Prony 
Modificado. No método de Prony, estimam-se as freqiiências 
através da solução de um sistema linear sobredeterminado e 
mais a solução de uma equação polinomial. Já no método de 
Prony modificado, as frequências são diretamente estimadas 
por meio de um algoritmo de otimização de Newton 


II. METODOLOGIA 


Neste trabalho utilizou-se um sistema de aquisição com- 
posto por um módulo condicionador de sinais e por uma 
placa de aquisição que é fixada internamente em um com- 
putador pessoal do tipo laptop. A placa de aquisição possui 
8 canais. A largura da banda de passagem do sistema é 
de 22,5 kHz. A frequência de amostragem pode ser fixada 
em até 500/c kHz, em que c é o número de canais. As 
amostras são digitalizadas em palavras de 12 bits. As medições 
foram realizadas conforme o esquema mostrado na Figura 
1. Os sensores de corrente funcionam por efeito hall como 
transdutores de corrente para tensão. A relação de transdução é 
de 100 mV/A. Os atenuadores de tensão possuem uma relação 
de transformação de 100 para 1. O condicionador de sinais é 
constituído por buffers e amplicadores que proporcionam um 
nível de tensão adequado para a placa de aquisição. 

O sistema de aquisição da National Instruments” M também 
inclui o programa Labview versão 7. O programa permite uma 
aquisição temporizada dos sinais da rede 


A. Armazenamento dos dados adquiridos 


A caracterização de harmônicos e inter-harmônicos no 
sistema elétrico deve considerar que existe uma incerteza 
intrínseca em relação ao valor do período fundamental do sinal 
analisado. Além disso, os componentes senoidais presentes no 
sinal são não estacionários devido ao comportamento dinâmico 
das cargas conectadas ao sistema. Obviamente em diferentes 
horas do dia ou em diferentes dias da semana a demanda de 
carga no sistema pode ter valores completamente diferentes. 
Em [13] sugere que é razoável um tempo de uma semana 
para caracterização de uma carga industrial em relação aos 
harmônicos ou inter-harmônicos. Contudo, se os sinais fossem 
adquiridos de modo ininterrupto durante uma semana, uma 
enorme massa de dados teria que ser armazenada na memória 
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Fig. 1. Esquema de medio dos sinais de tenso e corrente da rede eltrica. 


do laptop, o que inviabilizaria a análise dos dados. Assim, 
devido as restrições práticas, decidiu-se que janelas do sinal 
com duração de 1 segundo seriam adquiridas a intervalos 
regulares de 10 minutos. Assim, ao final de 1 semana de 
medições, seriam necessários 400 MB de memória para ar- 
mazenar as três tensões e três correntes de fase oriundas dos 
transformadores de corrente e de potência. Essa quantidade 
de memória é pequena em relação aos 20 GB de memória 
do disco rígido do laptop utilizado nas medições. Ao mesmo 
tempo, o intervalo de 10 minutos entre uma janela e outra é 
adequado pois o nível de distorção harmônica dos sinais não 
costumam sofrer alterações tão rápidas. Ou seja, é possível 
acompanhar a evolução nos níveis de distorções no decorrer 
da semana sem uma perda excessiva de informação. 

O intervalo de amostragem utilizado nas aquisições dos 
sinais foi fixado em 0,2 ms. Portanto, componentes senoidais 
de frequências de até 2500 Hz podem ser adquiridos corre- 
tamente, sem problemas de aliasing. De acordo com [13], 
componentes com fregiiências maiores são associados a am- 
plitudes muito baixas. Assim, mesmo que negligenciados, não 
afetam significativamente a estimativa de componentes com 
fregiiências menores. 


B. Processamento dos sinais 


O processamento dos sinais proposto neste trabalho com- 
bina o estimador linear de Kalman com o método de PM 
em uma análise off-line dos dados. Aqui, é importante salien- 
tar que devido a não linearidade do problema, a estimação 
das frequências pelo método PM é consideravelmente mais 
onerosa em termos computacionais do que a estimação das 
amplitudes pelo filtro de Kalman. Por isso, resolveu-se adotar 
um esquema de análise em que o método PM é aplicado 
em um número restrito de janelas obtidas durante todo o 
processo de aquisição. De maneira simplificada, a análise dos 
sinais pode ser melhor entendida com o auxílio do diagrama 
de blocos mostrado na Figura 2. O bloco A do diagrama é 
detalhado pelo fluxograma mostrado na Figura 3 e o bloco B 
detalhado pelo fluxograma mostrado na Figura 4. 

No fluxograma da Figura 3, ilustra-se a maneira pela qual 
o método de PM é aplicado aos arquivos disponíveis para a 
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Fig. 2. Diagrama de blocos simplificado da análise dos sinais adquiridos da 
rede elétrica. 


análise. Percebe-se que antes da aplicação do filtro de Kalman, 
o método de Prony modificado é aplicado em um conjunto 
restrito de arquivos, porém que abrangem todo o tempo no qual 
o processo de aquisição foi realizado. Assim, o método PM 
é sucessivamente aplicado aos arquivos cujas aquisições estão 
separadas por um tempo de 100 minutos. Devido aos arquivos 
serem nomeados segiiencialmente, o contador CONT deve 
ser incrementado de 10 unidades a cada laço do programa. 
Antes do método PM ser aplicado em cada janela, deve-se 
selecionar, entre as três tensões e três correntes de fase, qual 
sinal será analisado. As fregiiências estimadas pelo método 
PM devem ser armazenadas em um vetor vf, denominado 
vetor de fregiiências, a cada laço do algoritmo. Se em uma 
das iterações, o método PM fornecer uma ou mais frequências 
de valores já armazenados no vetor de fregiiências, então o 
algoritmo deve efetuar as médias dos valores próximos. As 
médias e os valores distintos obtidos obtidos pelo método 
PM, devem formar, juntamente com os valores previamente 
armazenados, um novo vetor de fregiiências. Assim se, por 
exemplo, na iteração i, vf; for expresso por: 


vfi = [60 120 180], 


e o método PM fornecer as frequências 30 e 60,5 rad/s, então 
na iteração à + 1, vetor vf;,1 será expresso por: 


vfiva = [30 60,25 120 180]. 


Observe que 60,25 é a média entre 60 e 60,5. Como saída, 
o fluxograma da Figura 3 fornece um vetor de freqiiência 
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Fig. 3. Fluxograma detalhando a aplicação do método de Prony modificado. 


que contem todas as fregiiências estimadas para diferentes 
janelas do grandeza selecionada (correntes ou tensões). Apesar 
das janelas estarem espaçadas de 100 minutos no tempo, a 
informação contida pelo vetor vf é significativa, pois o tempo 
total de aquisição da grandeza é bastante longo. 


O primeiro processo representado no fluxograma da Figura 
4 é necessário para evitar a polarização da estimativa das 
amplitudes pelo filtro de Kalman. A polarização ocorrerá se 
componentes senoidais de baixas amplitudes forem descarta- 
dos devido ao nível de thresholding estipulado pelo método 
PM. Para evitá-la e necessária a formação do novo vetor de 
fregiências vfn. O vetor vfn é a combinação de um vetor vp, 
de fregiiências previamente fixadas, e do vetor vf, oriundo 
das estimativas do método PM. Os elementos do vetor vp 
representam as frequências de 60 rad/s e suas 40 primeiras 
harmônicas. O vetor vfn é uma versão ampliada de vp. Assim, 
no processo de formação de vfn, incorpora-se ao vetor vp, 
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Fig. 4. Fluxograma detalhando a aplicação do filtro de Kalman. 


os elementos de vf que são distintos dos elementos de vp. 
Já os elementos de vp semelhantes aos elementos de vf, são 
substituídos pelos últimos. Assim, por exemplo, se o vetor vf 
for representado por: 


vf = [60,1 70 120,1], 
então o vetor vfn é expresso por: 
vfn = [60,1 70 120,1 180 240 -.. 2400], 


É claro que o conceito de semelhante ou distinto é rela- 
tivo e arbitrário. Na metodologia aqui proposta, considera-se 
frequências semelhantes aquelas cujos módulos das diferenças 
entre elas são inferiores a 0,5 rad/s. O primeiro processo 
minimiza a polarização na medida em que, pelo menos 
frequências harmônicas de alta frequência, estarão, necessaria- 
mente, consideradas nos regressores do estimador de Kalman, 
mesmo que sejam previamente descartadas devido ao nível 


de thresholding do método PM. É importante ressaltar que 
essas freqüências mesmo, não presentes no sinal analisado, 
não afetam a estimativa das amplitudes relacionadas com 
fregiiências, verdadeiramente, presentes no sinal. 

O segundo processo relacionado ao fluxograma da Figura 
4 é a aplicação do filtro de Kalman nas janelas do sinal 
analisado. Como já mencionado anteriormente, devido à lin- 
earidade, o problema da estimação das amplitudes é bem 
menos oneroso em termos computacionais do que o método 
PM, e assim o estimador de Kalman pode ser aplicado a todas 
as janelas adquiridas do sinal analisado. Isso significa que o 
contador CONT deve ser incrementado de uma unidade a 
cada laço principal do algoritmo. A aplicação do filtro de 
Kalman para cada janela do sinal resulta na evolução das 
amplitudes de cada componente senoidal presente em cada 
uma delas. As evoluções são denominadas pelos vetores: 


ec1, €ca, *::, €Cp, 


em que p é o tamanho do vetor vfn e ec; é o vetor que 
armazena a evolução do componente fundamental. Para se 
estimar a Distorção Harmonica Total (DHT) média relacionada 
a cada janela, é necessário, antes, obter-se a média dos vetores 
ec1, ecg, » €Cp. As médias são denominadas pelos 
escalares: mcy, mca, *-*, MCp. Dessa forma, o DHT médio 
para cada janela pode ser expresso por: 





ym + me3 +--+ me? 
DHT = 





(1) 


mc 


O DHT pode ser armazenado em um vetor, denominado 
DHT, a cada iteração do laço principal do algoritmo. Do 
mesmo modo, as médias de cada evolução podem ser ar- 
mazenadas nos vetores mc, mca, ---, MCp. Esses vetores 
podem ser salvos em arquivos e plotados para visualização e 
análise gráfica. 


II. RESULTADOS 


No primeiro exemplo mostrado neste trabalho, os sinais 
foram adquiridos na subestação de entrada de uma unidade 
industrial produtora de calçados com potência nominal de 10 
MVA. Na unidade estão instalados motores e conversores de 
potência. O esquema de como foram realizadas as medições 
está mostrado na Figura 5. Os sinais foram adquiridos no 
lado da alta tensão de um transformador de 13,8 kV para 
440 V através de três transformadores de potencial e três 
transformadores de corrente. A relação do transformador de 
potencial é de 120 para 1, e a relação do transformador de 
corrente é de 200 para 5. A aquisição dos sinais abrangeu 
um período de pouco mais de 5 dias. Na Figura 6, mostra-se 
o gráfico da evolução da Distorção Harmonica Total para a 
tensão da fase A durante o período de aquisição. Observa-se 
que a distorção não excede o limite recomendado pela norma 
IEEE 519 que é de 5%. Nas Figuras 7 e 8 mostram-se a 
evoluções da amplitude do quinto harmônico de tensão na fase 
A e do DHT da corrente na fase A, respectivamente. 

No segundo exemplo, os sinais foram adquiridos em uma 
subestação da Saelpa com uma linha expressa alimentando 
uma unidade industrial produtora de cimento. A aquisição dos 





Subestação de entrada - Planta Industrial 1 


TC 
(0) e 
q 


Trafo 
440/13,8 kV 
10 MVA A + 
Medição 
em painel 








sensor de corrente 











TP 
13.8 kV/115 V, 
S 


Y 


M 
Condionador 
de sinais 
Medição (y Atenuador 
em painel 























Y 
Placa de 
aquisição 























Fig. 5. Esquema de medição para a Planta Industrial 1. 
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Fig. 6. Evolução da DHT da tensão na fase A - Planta Industrial 1. 


sinais abrangeu um período de 8 dias. Na Figura 9, mostra-se 
o gráfico da evolução da Distorção Harmonica Total para a 
tensão da fase A, durante o período de aquisição. Na Figura 
10, mostra-se a evolução da amplitude do quinto harmônico 
de tensão na fase A. Observa-se que por volta do sexto dia de 
aquisições tanto o DHT quanto o quinto harmônico da tensão 
na fase A sofrem um aumento abrupto. A Saelpa reportou que 
houve um desligamento de energia neste dia. 

Também pode ser observado que no primeiro dia de 
medições o DHT de tensão na fase A ultrapassa o limite de 
5% recomendado pela norma. 


IV. CONCLUSES 


Neste trabalho foi apresentada uma nova metodologia de 
análise de harmônicos e inter-harmônicos de sinais de tensão 
e corrente na rede elétrica. A metologia supera as limitações 
da transformada rápida de Fourier e permite que a Distorção 
Harmônica Total dos sinais seja analisada por um longo 
período. 
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Fig. 7. Evolução da amplitude do quinto harmônico da tensão na fase A - 
Planta Industrial 1. 
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